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ABSTRACT 

A design study has been conducted for a direct air density measuring 
device that can be carried onboard an orbiting spacecraft in the altitude 
range from 140 to 280 km. The resulting design calls for a free molecular 
cryopumped gas collector of nearly unity capture coefficient, using an 
orifice collector and a cooled piezoelectric crystal acting as a micro- 
balance. Refrigeration is by heat sinks pre-conditioned before launch. 
Measurement is through beat-frequency (between cryopunq>ing crystal and 
a variable oscillator) converted to DC voltage, thence to telemetry. 
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TECHNICAL MEMORANDUM X-53468 


AN ORBITING DENSITY MEASURING INSTRUMENT 
SUMMARY 


A design study has been conducted for a direct air density measuring 
device that can be carried onboard an orbitirig spacecraft in the altitude 
range from 140 to 280 km. The resulting design calls for a free molecular 
cryopumped gas collector of nearly unity capture coefficient, using an 
orifice collector and a cooled piezoelectric crystal acting as a micro- 
balance. Refrigeration is by heat sinks pre-conditioned before launch. 
Measurement is through beat- frequency (between cryopumping crystal and 
a variable oscillator) converted to DC voltage, thence to telemetry. 


I. INTRODUCTION 


The need for the accurate prediction of spacecraft orbit decay has 
generated a need for an accurate determination of the pertinent atmos- 
pheric property - mass density, i.e., mass per unit volume, at the 
orbital altitudes of concern. The most direct way of measuring this 
quantity is to proceed by its definition, namely, to measure the mass 
of gas, regardless of composition, contained in a known ambient volume 
at orbital altitudes. Any technique attempting to make such a determina- 
tion must be capable of (1) collecting the mass of gas, (2) sensing the 
collected mass, (3) determining the volume from which the mass was col- 
lected, and (4) determining spatial position. 

Tracking provides the spatial position data, as well as spacecraft 
velocity; thus, an onboard instrument must consist of con^jonents capable 
of collection, sensing, and volume determination. 

Celestial Research Corporation has conducted a feasibility and design 
study to develop an air density measuring instrument based primarily on 
existing technology, so that the successful conclusion of the study would 
lead expeditiously to appropriate flight hardware without extensive 
research and development efforts. This report documents the results of 
the study on two possible methods and recommends one for hardware develop- 
ment. 


The nominal altitude considered is 180 kilometers; the vehicle 
assumed is one similar to the Saturn S-IVB. 
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II. BASIC SCHEME OF DIRECT DENSITOMETRY BY CRYOPUMPING 


The design approach is based upon the fact that a gas collection 
element aligned parallel to the spacecraft flight direction will sweep 
out an ambient gas volume 


V = A ut, (1) 

a o * ^ 

where Aq is the projected inlet area of the collector, u is the space- 
craft velocity, and t is the time of collection. 

The collection time is primarily a function of the mass sensing 
technique which in turn is determined by available instrumentation. Two 
approaches were investigated which more or less bracketed the possibil- 
ities, in that one instrument, a pressure transducer, relied on an inter- 
mittent collection cycle with pressure amplification, while the other 
instrument, a cryogenically cooled quartz crystal micro-balance, could 
continuously collect and sense the incoming mass . flux. Each of these 
measurement techniques had its advantages and shortcomings which will be 
discussed in detail in a later section. 

The scheme for the intermittent measurement technique involves the * 

collection of ambient gas for some specified time, the valving off of the 
collector with the captured cryogenically pumped gas, and the vaporiza- 
tion of the collected solid condensate by warming the collector to a 

known temperature. We can then determine the ambient mass density, , p^, « 

knowing the collector gas pressure as measured by the pressure transducer, 
the collector gas temperature, which is the same as the collector wall, 
the gas constant, %i, the collector measurement volume, Vm> and the ambient 
collected volume, from equation (1), 




Va Km ^ 


or 


V + V 

_£ a 

AoUt 
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where 


m = combined measurement volume 
c = collector 
g = pressure gage 
a = ambient. 

Pressure can be anq>lified to readily measurable levels by adjustment of 
the inlet area, the collection time, and the measurement volume. The 
measurement cycle is conq>leted by opening the collector inlet valve, 
thus venting the collected gas out to the free stream atmosphere and 
reducing the residual collector pressure to a negligible level relative 
to the pressure level during measurement. 

The continuous mass measurement technique is more direct in that a 
cryogenically cooled quartz crystal cryopumps the entering mass onto the 
resonating crystal which responds to changes in system mass by a shift 
in the resonant frequency; thus, mass flux is measured directly without 
intermediate measurements or calculations. The frequency shift is 
approximately a linear function of the collected mass until the shift 
becomes approximately 1 percent of the crystal's resonant frequency. 

For a crystal in thickness shear resonance, the sensitivity, i.e., the 
frequency shift with mass change, can be expressed as a function of the 
resonant frequency over the applicable linear range: 

Af = C Zto f|, (3) 


where 


C ~ constant 

= basic resonant frequency 
^ = collected mass 
AE = frequency shift. 

The question of whether or not the crystal need be thermally cycled to 
clean the crystal of accumulated condensate and thus avoid loading the 
crystal beyond its linear response range is influenced not only by orbital 
parameters but also by the required collector geometry for high capture 
coefficients and by the frequency shift sensing equipment. This will be 
discussed in detail in a later section. 



The refrigeration source and coolant system for the cryopumping of 
the collected mass is dependent upon the measurement technique. In the 
case of intermittent collection and pressure amplification, the entire 
collector wall is held at a low temperature during the collection cycle 
such that condensation of the entering gas molecules occurs with resultant 
condensate vapor pressure at least an order of magnitude below the ambient 
pressure. 

For air constituents, this means temperatures on the order of 20°K 
which can be obtained with either liquid hydrogen or helium as the 
refrigeration source. Once the gas sample has been collected, the col- 
lector wall, which has previously been functioning as a cryopump, must 
be raised to a temperature level which will again vaporize the gas sample 
for the pressure measurement. Once the pressure measurement is made, the 
gas sample is vented, and the collector is again cooled to approximately 
20®K and the inlet valve reopened for the next gas sample. This thermal 
cycling calls for a rather sophisticated refrigeration system which will 
be discussed later. 

In the case of continuous collection, the quartz crystal, which 
becomes the cryopumping surface, must be cooled to temperatures below 
20®K to collect the incoming mass. The collector walls are held at an 
intermediate temperature - high enough to avoid cryopumping yet low 
enough to reduce the energy level of impinging molecules to achieve 
eventual condensation on the crystal. This results in a relatively 
simple coolant system, since the heat load to both the collector walls 
and the crystal remains constant throughout the flight. The use of the 
crystal in this manner, i.e,, as a continuous mass collector, is con- 
tingent upon the ability of the sensing instrumentation to measure per- 
centage-wise very small frequency shifts, thus avoiding a thermal cycling 
of the crystal. 

The individual densitometer components will be studied separately 
in the following sections, and an overall system error analysis will be 
presented for the recommended system design. 


III. COLLECTION PROCESS 

A. Molecular Flux Measurement and Capture Efficiency 

The ti'-st step in this type of density meter is the collection 
phase. In this phase the free stream molecules are captured by condensa- 
tion on a cryopump. As a simple example, consider a cryopump consisting 
of a flat panel that is extended normal to the free stream (see Figure 1). 
If the flow is free molecular with respect to the cryopump, and if the 
molecules that impinge on the front face of the cryopanel are condensed, 
the mass will accumulate at a rate given by 



( 4 ) 
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• where 

N = molecular flux (molecules /sec-cm^) 

^ n = density (molecules /cm^) 

k = Boltzmann constant 
T = temperature 
M = molecular weight 
S = speed ratio = ■ ■ 

X(S) = exp (-S^) + n/jT S(1 + erf S). 

Uius, by measuring the mass that has been collected over a known period 
of time (t), an average free stream molecular flux can be determined, 

• If s, T and M are constant over that time period, the average free 
stream density can be determined from equation (4). 

If all of the impinging molecules are not condensed, equation 
V (4) must be modified. Assuming that only a fraction (q>) of the imping- 

ing molecules are captured, equation (4) becomes 


Nc =<pN, 


( 5 ) 


where 


cp = captured coefficient. 

Thus, in order to relate the measured collection rate and the 
free stream density, the capture coefficient must be known. If there is 
no interaction between the incoming molecules and the molecules that 
rebounded from the plate without being captured, the capture coefficient 
(qi) is equal to the surface sticking coefficient (o). The sticking coef- 
ficient is a more fundamental parameter since it is the probability that 
an incident molecule will be captured during a collision with a cryopun^>. 
In this application, the interaction between rebounding and incident 
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molecules is negligible (since the mean free paths involved are always 
much larger than the collector dimensions) so that in the case of the 
flat plate, <p = a. 

B, Surface Sticking Coefficient 

It is apparent from equations (4) and (5) that, for the simple 
flat plate configuration, the uncertainty in the measurement of the free 
stream density will be no less than the uncertainty in the sticking coef- 
ficient. The sticking coefficient is closely related to the critical 
energy level (Ec). This is the average maximum energy level that an 
incident molecule can have and still be captured during a collision with 
the wall. It is shown in Appendix A that there is considerable uncer- 
tainty in the value of the critical energy level, but it appears to be 
in the range of 1 to 25 times the binding energy of the molecule. The 
binding energy (Eb) is essentially the energy that must be expended to 
remove a condensed molecule from the condensate layer. This is comparable 
to the heat of vaporization or about .07 eV for a nitrogen molecule. In 
contrast, the energy associated with incident molecules (E^) during 
orbital conditions is about 9 eV. Since the incident energy is 4 or 5 
times the estimated upper limit of the critical energy level, it is 
doubtful that any significant number of the incident molecules will stick 
on this simple flat-plate cryopump. 

The negligible value for the sticking coefficient of incident 
molecules with orbital energy levels represents a fundamental limitation 
on the flat-plate cryopimip when operated in free molecular flow conditions. 

C. Capture Coefficient 

While the sticking coefficient is independent of the cryopump 
configuration, the capture coefficient is not, and can be increased. 

For example, consider a cryopump that consists of the interior of a 
cylinder and an end disc as shown in Figure 2. By aligning the cylinder 
with the free stream, the directed velocity of the molecules will cause 
most of the molecules to have their first collision deep inside the tube. 
While these molecules will probably not stick on this first collision, 
they will have a high probability of having several more collisions with 
the walls before passing back out the entrance. During each collision, 
the energy level will decrease, so that, if there are enough collisions, 
the energy level will fall below the critical level and the molecule will 
be captured. This is the same technique discussed by Wallace and Rogers 
[1] as a method for maintaining a low background pressure while test 
firing a rocket motor. 




The feasibility of this technique depends upon the energy lost 
per collision, or the accomnodation coefficient (a), which is the ratio 
of the energy lost in a collision to the maximum energy that could be 
lost. 


a = 


E. 

1 


- E 


off 


E. 

1 


F * 

wall 


( 6 ) 


Since there is a lack of e:q>erimental data for the accommodation 
coefficient of a high energy gas iiiq)inging on a condensate layer, it is 
necessary to extrapolate the available data using theoretical results as 
a guide. As outlined in Appendix B, this approach yields an estimate of 
CC = .77. This is asstimed to be independent of the molecular energy level, 
so that the number of collisions (v) required to lower the energy level 
to the critical value can be obtained from 


(1 - a)"' 



( 7 ) 


When this is evaluated with a = .77, = .07 and Ei = 9, it is found 

that V = 3.3. Thus, 3 or 4 collisions would be required to lower the 
energy level to the point that the molecule would stick on the next col- 
lision with the cold wall. 

D. Cylindrical Collector 

An analysis based on a Monte Carlo program developed by Ballance 
[2] has been made of a series of cylindrical cryopumps operating at S = <» 
and S - 10. In using this approach, it is assumed that all molecular 
rebounds are diffuse. Typical results are shown in Figure 3. It is seen 
that, when S = “, there is a continuing benefit from increasing the L/D 
of the cylinder. This is to be expected, since for S = oo, all the incident 
molecules strike the aft disc; therefore, increasing the tube length 
increases the distances the molecules must travel to the opening. Increas- 
ing this distance increases the probability that a molecule will have 
enough collisions to decrease the energy below the critical energy level 
so that it will be captured. 

Increasing the L/D of the cylinder is also effective when S = 10 
if L/D is small. However, if the L/D is large, there is no significant 
benefit from further increases. In the case of a large L/D, the capture 
coefficient is limited by the loss of those molecules that initially strike 
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the cylindrical walls near the inlet. Since this loss is essentially 
independent of the L/D of the cylinder, there is no benefit in increas- 
ing the ratio beyond L/D = 6. From Figure 3 it appears that there is an 
uncertainty of at least 5-10 percent in the value of the capture coef- 
ficient for a cylindrical cryopump. 

Consideration of the losses due to sputtering of the condensate 
layer by the high energy incident molecules also increases this uncer- 
tainty. No experimental data are available for the sputtering of an air 
condensate layer by incident air molecules, but a crude correlation (see 
Appendix C) based on results of metallic sputtering indicates that the 
yield would be of the order of unity; i.e., each high energy incoming 
molecule would sputter one condensate molecule. While it appears that 
neither of these molecules would have sufficient energy to cause further 
sputtering, this initial sputtering increases the uncertainty of the 
capture coefficient by about 5 percent. 

Since the estimated uncertainty is quite high and becomes even 
higher if angle of attack effects are considered, this cylindrical con- 
figuration does not appear satisfactory for use in an orbital collector, 
unless the uncertainties can be further reduced by some extensive experi- 
mental investigation. Since there are no facilities that can properly 
duplicate orbital conditions, it is advisable to find a configuration 
that is less sensitive to the uncertainties in sticking coefficient, 
accommodation coefficient, and sputtering yield. 

E. Orifice Collector 

The capture coefficient based on the inlet area can be further 
increased by using an inlet diameter that is smaller than the tube diam- 
eter (Figure 4), For example, if the inlet diameter is 1/3 of the tube 
diameter, the molecular loss rate is reduced to about 15 percent of that 
obtained when the inlet diameter is equal to the tube diameter. This 
approach has the further advantage that the collector is much less sensi- 
tive to angle-of-attack effects. On the other hand, the reduced inlet 
area results in a lower flow rate into the collector, so that for a fixed 
altitude either the collection time or the gage sensitivity must be 
increased. These factors will be considered in more detail in another 
section. 


While this reduced inlet area configuration does decrease the 
uncertainties in the capture coefficient, the configuration shown in 
Figure 5 reduces the uncertainties even further. In this configuration, 
the molecules that enter through A© impinge on the wall As that is cold 
but not cold enough to cryopump. After a sufficient number of collisions 
with this cold wall (T < SO^'K), the energy level of the molecules is 
below the critical value and the molecules will condense upon their next 
collision with Ac, the cryopumping surface (T < 20®K), During the 
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collisions that are required to cool the molecules below the critical 
energy level, some molecules will pass out through the inlet. If the 
rebounds are diffuse, the fraction of the molecules lost per wall col- 
lision is Aq/As* If there are v collisions required to cool the mole- 
cules, the total losses during this phase will be vAq/As. 

There will also be molecular losses during the condensation 
phase. These fractional losses will be approximately equal to Aq/Aq. 

Thus the total inlet losses during collection will bevA^/Ag + A^/A^,. 

To minimize errors, it is necessary that these losses should not exceed 
1 to 2 percent. This requires that Aq/Aq and vAq/A^ each be about 0,01, 

The value of v depends upon the thermal accommodation coefficient. 
This in turn varies with the intermolecular force between the gas and the 
wall, and the ratio of the molecular weights. If the system is fabricated 
from aluminum, the wall molecular weight is comparable to that of the col- 
lected gas. Furthermore, the intermolecular force between the gas and the 
metal wall should exceed the corresponding forces between the gas and a 
condensate layer; thus, the gas-wall accommodation coefficient should 
exceed the gas -condensate accommodation coefficient. Therefore, it 
should be conservative to use the value for the accommodation coefficient 
obtained previously ((X = 0.77). It then requires 4 collisions to reduce 
the molecular energy below the critical value. This then requires that 
4Ao/As <0.01. 

In general, since this condition can always be met, the losses 
will be dominated by the ratio Aq/Aq. 

A further advantage of this approach is the reduced possibility 
of sputtering of the condensate layers. This is due to the initial col- 
lision for essentially all the incoming molecules being with a cold metal 
wall and not a condensate layer. The possibility for sputtering on sub- 
sequent collisions is also small because of the low average energy level 
of the molecules striking the condensate layer. 

While the basic analysis has been made for a spherical configura- 
tion, it may be more convenient to use a cylindrical configuration. This 
change would not alter the results, since in general the inlet area is too 
small to maintain the pressure gradients that are required to prevent a 
nearly uniform distribution from being established. 

F, Noncondensable Gases 

A fraction of the ambient gas sample captured by the collector 
will be made up of gas species which are either noncondensable at the 
20®K cryopump temperature or will have an unacceptably high solid vapor 
pressure. This fractional amount will constitute an inaccuracy in the 
determination of the ambient density by the proposed mass collection 





principle. The three gases of concern, helium, hydrogen and neon, are 
present in small amounts in the homosphere where the various species are 
distributed according to a common number density-altitude profile due to 
general turbulent mixing. Above the mictopause, however, where the homo- 
sphere ends and the heterosphere starts, changes in the atmospheric com- 
position occur as a result of gravitational separation and by interaction 
with incoming radiation. The gravitational separation tends to increase 
the relative percentage of the light gases in the atmosphere and thus 
potentially aggravates the noncondensable gas problem. Unfortunately, 
almost no data are available concerning the composition in the subject 
range of interest (140 to 280 km altitude); however, theoretical esti- 
mates can be used to arrive at approximate effects. If 100 km is con- 
servatively used as the mictopause and the composition is calculated for 
the extreme case of 280 km altitude, helium is found to constitute 0.014 
percent of the entering gas sample mass flux and hydrogen is 0.002 per- 
cent. Neon, being a heavier gas, varies by only a small amount from the 
homosphere percentage, i.e., 0.0012 percent. The presence of noncondens- 
able gas species in the entering gas sample thus constitutes an extremely 
small error in the ambient density determination by mass collection. 


IV. MEASUREMENT TECHNIQUE AND ACCURACY 


A. Intermittent Collection 

1. Minimum Collection Time 

The minimum collection time may be set by the desired 
accuracy in measuring density variations as a function of altitude or 
by anticipated diurnal density variations. In either case, the optimum 
configuration of collector and pressure transducer should produce the 
maximum measurement pressure for a given collection time. From equation 
( 2 ), 


P = 
m 


A ut R 
o m 

^a V + V 


m 


This can be expressed in terms of the important collector 
parameters as follows; 


P 


m 


Pa ^m ^m (VDc)^ 


O), 


( 8 ) 
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where 


(V^./Vg)2/3 

P " (1 + V^/Vg) • 


* The collector length-to-diameter ratio (L^/D^) and the inlet 

orifice-to-collector diameter ratio (Dq/D^.) are determined by requirements 
of capture coefficient; thus, the measuring pressure (P^) is maximized by 
the function 3. Maximizing 3, we find that a collector-to-gage volume 
ratio of 2 will yield the maximum pressure, although volume ratios between 
1 and 4 give results within 95 percent of the maximum. 

This optimum pressure is plotted versus the collection time 
normalized by the orbital period in Figure 6. In this case, the collec- 
tor- tube- to-aperture area ratio is taken as 10 with the collector L/D =5. 
This geometry gives an estimated capture coefficient of 98 percent as dis- 
cussed in the previous section. The most promising pressure transducer 
encountered in the study (MKS Baratron gage) has a full range pressure 
measurement capability of 1000 microns and a volume of 18 cm^. An evalua- 
tion of the minimum gage measurement pressure thus relates the minimum 
collection time and the measurement altitude. 

% 

To resolve the diurnal density variation, the collection 
time should not be greater than one- tenth of an orbit, i.e., t/t^ < 0.1. 
Available data indicate that the magnitude of this diurnal variation may 
^ be on the order of a factor of two in measured density. This represents 

a significantly greater density change than that encountered due to an 
orbital decay. Approximately 40 orbits would be required, for instance, 
to obtain this variation from strictly altitude change in the case of the 
S-IVB stage of the Saturn vehicle decaying frcm an initial 180 km circular 
orbit. Thus, the controlling time requirement is the diurnal density 
variation; this limits the collection time to not more than 0,1 orbit. 

2. Altitude Range 

The maximum altitude at which meaningful density determina- 
tions can be made is related to the minimum measurement pressure. If 
this minimum pressure is assumed to be 10 percent of the instrument’s 
full range, i.e,, approximately 100 microns. Figure 6 indicates a maxi- 
mum measurement altitude on the order of 220 km due to insufficient 
pressure. 
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The low altitude limit is set by gas dynamic considerations 
rather than the measuring gage, since the assumption of free molecule 
flow at the collector inlet is basic to the determination of the swept 
ambient gas volume. The free molecule flow limitation is characterized 
by the diameter of the vehicle and the mean free path for collisions 
between free stream molecules and those emitted from the vehicle. 
Scattering of the free stream molecules by the emitted molecules begins 
to become significant, that is, greater than 10 percent, when the ratio 
of diameter to mean free path approaches unity. The ratio of free stream 
mean free path to the mean free path between emitted and free stream 
molecules at orbital velocities is about 12. Therefore, for a vehicle 
diameter of 6.6 meters, the required free stream mean free path is about 
79 meters which corresponds to an altitude of 165 km. This lower limit 
can be reduced to 140 km by applying first-order corrections to the 
measurements to account for the scattering effects of the free stream. 
First- coll is ion scattering calculations of the type developed by Wain- 
wright and Rogers [3] and others can be made for mean free path lengths 
down to 1/4 of the 79 meter threshold. The design operational range of 
the instrument, then, is 220 km to 140 km, and spans a factor of approxi- 
mately 17 in density. 

3. Measurement Accuracy 

The overall accuracy of the measurement system will depend 
upon the confidence which can be placed in the primary measurements of 
pressure, temperature, vehicle speed and altitude and in the evaluation 
of certain systematic effects, such as the collector capture coefficient, 
distortion of the free stream flow field by the vehicle, leaking of the 
measurement chamber, etc. All these factors are represented in the equa- 
tion used to calculate the density. 


P V + V 

Q = I- . — s 

^a ^ R T uA tcp ’ 
mm o ^ 


(9) 


where is the ambient density of the atmosphere, Pjjj and T^ are the 
measured pressure and temperature in the collector, ^ is the gas con- 
stant of the gas in the collector, Kl is a collector leak factor, 

(Vg + Vg)/AQ is the ratio of collector and gage volume to collector 
aperture area, cp is the capture coefficient, u is the free stream velocity, 
and t is the collection time interval. 

The collector pressure and temperature are the primary 
measurements, and their accuracy will depend upon the quality of all 
elements in the measurement chain from the transducer to the ground- 
based telemetry station. Attention is directed in this study to system 



elements up to the telemetry encoder and transmitter with the view that 
the overall system accuracy up to this point will determine the quality 
of standard telemetry elements which can be justified* Mechanical dia- 
phragm type pressure transducers in the 1000 micron pressure range are 
presently available which can produce 0*5 percent measurements* Deflec- 
tion of the diaphragm is detected by variations in capacitance which is 
produced relative to a fixed surface* The transducers are sensitive to 
temperature gradients; but since the entire densitometer instrument will 
» be conditioned and held at liquid nitrogen temperature, transducer tempera- 

ture sensitivity should not pose a problem* The volume- to-aperture area 
ratio, (y^ + Vg)/AQ has an uncertainty of 3 percent. 

It is anticipated that the collector chamber temperature 
will be controlled rather than measured. Control would be achieved by 
regulating the pressure over boiling liquid nitrogen. At 100°K, a 7 per- 
cent tolerance in the regulated pressure yields only a 1 percent error in 
temperature. In addition to this advantage, boiling point control elimin- 
ates the need for one telemetry channel. 

Determination of Rj^, the gas constant of the collected gas 
at the measurement temperature, Tm, involves some estimate of atmospheric 
composition at the measurement altitude and an assumption that all decom- 
posed molecules collected in the atmosphere will be recombined by the time 
of measurement in the chamber. The latter assumption is well founded 
^ because of the high collision frequency of the captured gas at the measure- 

ment pressure relative to the time required to evaporate the collected gas 
and measure. 

^ The errors induced by uncertainty in composition depend upon 

the molecular weights of the components and the composition. If, for 
example, the concentrations of nitrogen and oxygen are considered with 
the concentration of oxygen equal to 0, the error in evaluating Rm rela- 
tive to the uncertainty in 9 is given by 


dR/R I , (10) 

e(Mo - %) 


where and Mq are the molecular weights of nitrogen and oxygen, respec- 
tively. It is seen that a 20 percent uncertainty in 0 yields only a 
0.55 percent error in Rm. Thus, a relatively crude estimate of oxygen- 
nitrogen composition will suffice to give a value for Rm whose accuracy 
is consistent with the other factors of the density determination. 
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The leak factor, Kj^, is associated with the integrity of 
the pressure seal for the measurement chamber enclosing the collector 
during the warmup and measurement cycle. The only seal in the gas col- 
lector system is in the shutter over the inlet aperture. For the standard 
flight conditions associated with a 180 km orbit and assuming warmup and 
measurement time of the same order as the collection time interval, an 
uncertainty of 8 x 10"^ liter/sec in leak rate would yield an uncertainty 
in measured pressure of 1 percent. Leak rates of at least an order of 
magnitude higher could be tolerated provided it is known and repeatable 
to within 8 x 10“® liter/sec. It is anticipated that a teflon 0-ring 
type seal can be developed to meet these requirements. Tracking pro- 
cedures are capable of providing measurements of vehicle velocity to 
within 0.1 percent. The shutter governing the collection time, t, would 
be controlled by a programmer which would govern all the functions of the 
instruments. It is estimated that the opening and closing time would be 
of the order of 0.5 sec with a repeatability of 0.05 sec. Thus, with a 
normal collection time of 500 sec, the accuracy of t should be 0.1 percent 
or better. 


The last major item which must be known for the evaluation of 
density is the capture coefficient, cp. If reference is made to Section 
III on the collector design, it is seen that the cylindrical tube collector 
equipped with the proper aperture is capable of capture coefficients on 
the order of 98 percent. Thus, the fraction of swept molecules which 
are lost, i.e., the complement of cp, is approximately 2 percent. 

Intermittent collection for the determination of ambient 
density by pressure transducer measurement thus has an uncertainty of 
approximately 8.2 percent up to the telemetry system which in itself has 
approximately a 1 percent uncertainty. 

B. Continuous Collection 

1. Maximum Collection Time 

A resonating quartz crystal responds to changes in system 
mass, e.g., the deposition of a condensing gas, by a shift in the resonant 
frequency. The crystal can therefore be employed as a micro-balance, if 
other potential indeterminables are eliminated. Certain crystal cuts are 
highly sensitive to temperature; however, cuts can be made which eliminate 
this effect over the desired temperature range. From the work of Waters 
and Raynor [4] and Stephens [5 and 6], a 39®40' A-T cut has been found to 
be insensitive to temperature changes below 20®K. The sensitivity of 
the crystals in thickness shear resonance as a function of resonant fre- 
quency was presented in equation (3). 



M. = CAmf^ 
r 

where 

C = 2.26 X 10*^ 

Af - cps 
Am - gms/ cm^ 

- cps. 

The crystal response becomes nonlinear when the frequency shift due to 
mass addition becomes approximately 1 percent of the resonant frequency. 
The allowable mass change for linear cyrstal response is thus 


Am 


max 


4.42 X 10^ 
f 


( 11 ) 


V 


Evaluating the system at 180 km altitude from available 
atmospheric property data then gives a maximum continuous collection 
time for linear response: 


where 


t 

max 


lO^Q 

(A /A )f 
^ o c^ r 


secs , 


Aq = aperture area 
Ac = crystal area. 


( 12 ) 


An earlier discussion indicated the necessity of area ratio 
values, (Ac/Aq), of approximately 100 to minimize measurement error due 
to escaping molecules; thus, the normalized maximum continuous collect- 
ing time becomes 


t 

max 



1.8 X 10^ 

f 

r 


(13) 
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where again to is one orbit period - in this case 5600 seconds. Crystals 
are available in the 1 to 10 megacycle range which are adequate for the 
orbit lives under consideration. The linear response range of operation 
is not an absolute limitation, of course, and calibration of the crystal 
response over a wider frequency shift range would extend the data-produc- 
ing collection time. The disadvantage of this scheme would be the incon- 
venience of summing the accumulated mass in order to determine the changing 
crystal sensitivity during the flight. 

2. Altitude Range 

The maximum altitude at which continuous collection for 
density measurement with a quartz crystal can be used is related to the 
minimum measurable condensed mass change and the length of time between 
measurements. The frequency shift, equation (3), can be related more 
directly to collector parameters and the collection time by rewritting 
as follows: 



(*/Aq) ^ 


o 


or 



o* 


(14) 


From earlier discussions certain restraints are imposed 
upon some of these parameters. It was established, for instance, that 
the area ratio, Aq/Ac» must be on the order of .01 for a high capture 
coefficient, and t/to can be no larger than 0.1 in order to measure 
diurnal density variations. The orbital period, to, is not greatly 
affected by orbit altitude nor is the vehicle velocity, u, which together 
with the ambient density, p^, forms a product which is the entering mass 
flux per unit area, m/Ao- It is apparent then that the minimum measurable 
value of ambient density, p^, is simply a function of the accuracy with 
which the frequency shift, AE, can be measured. The percentage frequency 
shift for a given mass change is proportional to the basic crystal resonant 
frequency; 


(15) 


f 


= CAm f 
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Therefore, selection of high crystal frequency increases the measuring 
sensitivity, although as shown in the last section this implies a shorter 
continuous collection time before nonlinear crystal response occurs. 

Since the frequency shift sensing and telemetry system has 
a limiting accuracy of approximately 1 percent, it is important that the 
total frequency shift be limited as much as possible in order to minimize 
inaccuracies in p . In the case of continuous collection, however, this 
total frequency shift may be on the order of 10,000 - 100,000 cps for a 
normal orbital lifetime. To reduce the magnitude of the shift, either 
the crystal must be thermally cycled to revaporize condensed mass between 
measuring points, thus returning the crystal to a zero frequency shift 
condition, or a beat frequency system must be employed in such a manner 
as to reduce the shift. The thermal cycle scheme introduces a complexity 
in the cooling system which is undesirable. On the other hand, the beat 
frequency technique increases the sophistication of the electronics sens- 
ing system, although, in this particular application, the cooling system 
is by far the more difficult component system. 

Smaller frequency shifts can be obtained by mixing the out- 
put of the measuring crystal with that of a variable reference oscillator 
The resulting beat frequency is kept within a prescribed frequency range 
by driving the reference oscillator to new frequencies as the beat fre- 
quency reaches maximum set values . The reference frequency is varied by 
means of a frequency-voltage converter operating a voltage-controlled 
drive motor on the variable oscillator. For the particular flight orbits 
considered in this study, a converter providing a 0,5 to 4.5 volt DC 
signal corresponding to a 50 cps to 450 cps beat frequency is appropriate 
This sensing system is shown schematically in Figure 7. 

The maximum altitude at which the proposed instrument can 
measure density thus becomes a question of the desired accuracy of the 
measurement. The telemetry system results in a 5 cps uncertainty for 
the beat frequency system; thus, a frequency shift of 50 cps yields 10 
percent uncertainty in the density determination. Equation (14) thus 
indicates the following minimum measurable densities as determined by 
the crystal resonant frequency and assuming the use of the aforementioned 
beat frequency system: 


V 

f 

r 



p 

h 



1 me 

5 

x 

10"® kg/m^ 

134 

km 


5 me 

2 

X 

kg/m^ 

220 

km 


10 me 

5 

X 

10-^^ kg/m^ 

283 

km 
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The low altitude measurement limitation can be imposed by 
the same gas dynamic condition which sets the limit for the intermittent 
system, or it may be determined by overloading of the crystal by the high 
mass flux at the lower altitudes. Based on the assumed flight vehicle, 
the lower altitude limit set by deviation from free-molecule flow condi- 
tions was found earlier to be 140 km. At this altitude, the maximum con- 
tinuous collection time is 


3.2 X 10^ 

tn f. 


Thus, a 5 me crystal would remain in the linear response 
regime for 6 orbits, and data could be obtained for continuing orbits. 

It must be concluded that mass flux is not the determining factor, and 
that the low altitude limit is 140 km, resulting from the vehicle size 
influence on the departure from the free molecular flow condition. 

3, Measurement Accuracy 

The overall accuracy of the continuous collection measure- 
ment system depends upon the factors represented in the equation used to 
calculate the ambient density by this technique: 


^ ^ 

~ Cf^ ucp(Ao/Ac) (t/t^) t^, 


(17) 


The error in density determination due to the error in fre- 
quency shift measurement is d therefore, the magnitude of the error 

is determined not only by the frequency sensing equipment accuracy, d Af , 
but also by the magnitude of the frequency shift itself. The uncertainty 
is thus influenced by the collection time for a particular measurement 
and by the orbit altitude. The frequency- to-voltage converter discussed 
previously is the only component in the best frequency sensing system up 
to the telemetry encoder which can introduce reading errors. This type 
of converter is a standard commercial component with 0.1 percent full 
scale accuracy. There will, however, be errors on the order of 1 percent v 

in the telemetry of the beat frequency generated voltage signal. 

Since the maximum beat frequency is a set value, the 
inaccuracy cannot be less than 1 percent and can be more for short 
measurement intervals or for higher altitudes which result in less than 
full scale frequency shifts. Evaluation of the frequency shift from 
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equation (17) with a 10 me crystal and an orbit altitude of 180 km indi- 
cates that the full scale shift of 500 cps will occur in approximately 
1/12 of an orbit which is less than the maximum allowable time. Thus, 
at 180 km which has been considered to be the nominal orbital altitude 
in this study, the frequency shift has an uncertainty of 1 percent. 

In this type of sensing system, the actual value of the 
reference oscillator frequency is immaterial, and it is only Important 
that the stability of the frequency during a particular measurement time 
period be assured. Temperature control of those standard oscillators 
results in extreme stability, and frequency shifts in the reference 
oscillator will be negligible. 

As previously discussed, tracking procedures provide a 
measurement of vehicle velocity to within 0.1 percent. The fraction of 
swept molecules lost due to non-unity capture coefficient is approximately 
2 percent. 


The inlet aperture to crystal area ratio, Aq/Ac, has an 
uncertainty of 3 percent. The crystal* s resonant frequency determination, 
as well as the collection time between measurements, contributes a negli- 
gible uncertainty. 

Continuous collection for the determination of ambient 
density by a cryogenically cooled quartz crystal micro-balance thus has 
an overall system measurement uncertainty of approximately 6.1 percent, 
including telemetry error. 


V. COOLANT SYSTEM 


A. Intermittent Collection 
1. Cooling Scheme 

The component elements of the intermittent collection system 
are shown schematically in Figure 8. The major elements are the collector, 
the pressure transducer, the aperture shutter, the two thermal switches, 
the heat transfer rod, a liquid nitrogen source, and the liquid helium 
reservoir. A typical measurement cycle begins with the cooldown of the 
collector to liquid helium tenq)erature by closure of the transfer rod 
thermal switch. During this cooldown, the aperture shutter is closed 
and the collector is isolated from the liquid nitrogen source. After 
collector cooldown has been achieved, the shutter will open for the pre- 
scribed collection time, permitting a known volume of ambient air to be 
captured. The next step is the vaporization of the condensed air sample 
and pressure measurement by the transducer. This is accomplished by 
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opening the transfer rod thermal switch and bringing the collector into 
contact with the liquid nitrogen source thereby warming the collector to 
100®K. The shutter is closed and sealed against the inlet aperture during 
the warmup to liquid nitrogen temperature. After warmup is achieved and 
the pressure measured, the shutter is opened, venting the collected gas 
out to the atmosphere and reducing the residual pressure to a negligible 
level. The collector is then ready for another liquid helium cooldown 
and sample collection. 

2. Refrigeration Requirements 

The quantity of liquid helium required to accomplish the 
measurement mission depends upon the requirements per measurement cycle 
and the total number of cycles in the mission. A frequency of 10 measure- 
ments per orbit was previously found to be required to resolve diurnal 
density variations. The descent from 180 km to 140 km will require 55 
orbits of the assumed flight vehicle. Therefore, a maximum of 550 measure- 
ments would undoubtedly suffice for both altitude variation measurement 
and the diurnal variation. Therefore, the number of measurement cycles 
could be considerably reduced if the operation of the aperture shutter is 
controllable. 


(a) Cooldown Requirement 

The liquid helium volume requirement for collector 
cooldown is related to the collector heat capacity and a specific liquid 
requirement as follows; 




5 P I 
c^c 


(V ) 

g 


2 . 1 3 


(Vc/Vg)^/^(Ac/Vc)^/=", 


(18) 


where 

6 = collector wall thickness 
= collector density 
= liquid helium density 

A = total collector area 
c 

I = specific liquid requirement. 


From this equation, the collector L/D for the minimum value of 
can be found which will lead to the minimum helium requirement. Solution 
in this manner results in an L/D = 1 collector as the optimum configura- 
^ tion. Previous discussions indicated, however, that an L/D = 5 cylindrical 

collector was required in order to achieve near-unity capture coefficients. 
Evaluating this configuration as cou^ared to the optimum configuration 
indicates only a 15 percent penalty in the helium requirement for an 
L/D = 5 collector. 

The specific liquid requirement, is the number of 
grams of liquid helium required to cool a gram of the collector material 
over a specified temperature interval, and has been analyzed by Jacobs 
[7]. In this particular application, copper would probably be used for 
the collector walls because of high thermal conductivity characteristics 
as low temperatures. A value of | = 0.5 applies to the cooling of 100°K 
copper to liquid helium temperature, if only the heat of vaporization is 
available. Use of a boiloff vent system which makes use of the gaseous 
helium's heat capacity in the cooldown would reduce | to 0.03. The more 
conservative figure will be. used in the following calculations. 

The optimum collector to gage volume ratio, Vc/Vg, was 
earlier found to be 2 and the gage volume is on the order of 18 cm3. For 
a collector wall thickness of 5c = 0.005 cm, the required liquid helium 
volume for cooldown is thus calculated to be 12.4 cm^ for each cycle. 

% 

(b) Condensed Gas Requirement 

During the measurement cycle, the collector walls and 
> transfer rod transfer the heat load associated with the condensation of 

the mass of entering free stream air to the liquid helium where boiloff 
will balance the heat load. The energy to be removed is essentially the 
kinetic energy of the entering air and can be related to the liquid helium 
volume requirement as follows: ^ 


where 




He 

%e 




H 


V 


Pa 


> 


= helium heat of vaporization. 


(19) 
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Using t/to = 0.1 for the sampling time and the optimum collector to gage 
volume ratio, 


^ = 3.71 X 10^2 p^. (20) 

He 


Evaluation at a nominal 180 km altitude results in a 
requirement of 2.22 cm® of liquid helium for each measurement cycle. 

This volume requirement is directly proportional to the ambient density 
and thus is easily evaluated as a function of altitude. 

(c) Helium Reservoir 

The requirement of 14.6 cm®/ cycle of liquid helium for 
a maximum of 550 cycles Implies a helium reservoir volume of 8 liters. 

In addition to the heat load associated with the gas collection process, 
there are heat leaks through the insulation of the reservoir. Assuming 
the use of multilayer insulation and 80®K surrounding temperature results 
in a requirement of 0,2 cm®/ cycle of liquid helium. To the total cycle 
requirements should be added an additional 16 cm® of liquid helium to 
account for reservoir insulation losses during a 12-hour prelaunch period. 
The total amount of helium required for 550 measurement cycles plus a 
12-hour hold period is 8166 cm® which implies a 10-inch diameter spherical 
reservoir. 


B. Continuous Collection 
1. Coolant Scheme 

Using the quartz crystal sensor as a continuous cryopump 
for the entering mass flux reduces the refrigeration requirements con- 
siderably compared to the intermittent collection system, inasmuch as 
there is no need for a thermal cycling phase in the coolant system, thus 
eliminating a major heat load. 

The reduced level of heat load makes possible the use of 
a precooled metal mass acting as a heat sink for a simple source of 
refrigeration.* The coolant system, schematically shown in Figure 9, 


An alternate cooling scheme by hydrogen sublimation is presented in 
Appendix D. 
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consists of the crystal heat sink mass (nominally at liquid helium tem- 
perature) in an insulated container, a conduction path from the heat sink 
to the crystal, and a gas collector constructed such that the walls act 
as a heat sink at a temperature near 100°K. The two heat sinks are pre- 
cooled before flight with liquid helium and liquid nitrogen by means of 
coolant coils integral to the heat sink masses. The collecting period is 
initiated at the desired altitude by permanently removing an aperture cap 
which has previously maintained vacuum conditions in the collector, thus 
e keeping the heat load to a minimum during prelaunch holds and the initial 

portion of the flight. The heat sinks slowly rise in temperature during 
the measurement period. Accurate measurement conditions are maintained, 
however, as long as the crystal temperature remains at a level which 
results in condensate vapor pressures at least an order of magnitude less 
than the gas sample pressure, and the collector walls remain below a 
level which would produce an unacceptable radiation heat load to the 
crystal or which would not reduce sufficiently the temperature of incom- 
ing gas molecules. The analyses in the following sections result in a 
crystal heat sink of 55 cm^ of lead and a collector wall heat sink of 
113 cm^ of aluminum. 

2. Refrigeration Requirements 

A hypothetical mission will be used in this analysis of the , 
required refrigeration system, consisting of an unattended 12-hour hold 
% period of the vehicle before launch followed by a 55-orbit flight. 

(a) Crystal Sensor 

> (1) Radiation Heat Load - The crystal receives a 

radiation heat load from the collector walls not only during the measure- 
ment period but also during any unattended prelaunch hold periods. The 
design of the crystal heat sink mass must take into account, then, the 
operational problems involved in large vehicle launches. 

The collector wall being a heat sink will be varying in 
temperature during the flight; however, a conservative heat load is obtained 
if a wall temperature of 90®K is used in the calculation. The crystal can 
conservatively be considered to have an absorptivity of unity and the col- 
lector wall will be polished to an emissivity of 0.1. 

The maximum available crystal diameter is approxi- 
mately 2.2 cm; therefore, A^ = 3.95 cm^, and from the radiation heat 
transfer equation. 


= (5.67 X 10-12) . j4) 


and 


( 21 ) 


Qj. = 1 milliwatt. 
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(2) Gas Condensation Heat Load - The entering mass 
flux is reduced in temperature by the cold collector walls before 
encountering the crystal. Assuming the gas temperature is reduced to 
100®K, the energy removal required, Zlh, for condensation on the crystal 
is 83 cal/gm. 


Using an altitude of 140 km for a conservative mass flux. 


hi/Aq = 26.5 X lO"”^ gm/sec - cm^. 


The condensation heat load, using the fact that 
Ac/Aq = 100, is then 

^cond ^ (*/Ao) Aq = 0.037 milliwatts. 


(3) Crystal Temperature Gradient - The equation for 
the temperature gradient in a flat conducting disc with a uniformly 
imposed heat load is 


T 


r 



- r^) 
4k6 


( 22 ) 


where ' 

q = heat load 

k = thermal conductivity of disc 
& = disc thickness 
R = edge radius 
r = local radius. 

The total imposed heat load due to radiation and gas > 

conduction is 


q = 2.6 X 10"^ milliwatts/cm^. 

Therefore, for a 0.015 cm thick quartz crystal 
(k = 5 watts/ cm- “K), the maximum temperature difference caused by the 
imposed heat load is AT = 0.005°K. 
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(4) Crystal Heat Sink - The material used in the heat 
sink should have a high specific heat in the applicable range of tempera- 
ture to minimize heat sink weight, as well as a high density which will 
tend to minimize the displaced onboard volume. Lead is an appropriate 
heat sink material on this basis, and the following equation can be used 
to calculate the required mass of metal required for a given time interval 


/ 


m = 


Q / dt 




(23) 


The allowable heat sink temperature rise is determined 
by the criterion that the condensate vapor pressure should be an order of 
magnitude less than the collector chamber pressure for negligible error 
in condensed mass determination. With an area ratio A^/Aq = 100, and with 
with a speed ratio of approximately S = 10, the crystal has approximately 
ten times the volumetric pumping speed of the inlet aperture. The col- 
lector pressure during measurement will therefore be approximately 1/10 
the ambient pressure. Using the condensate vapor pressure criterion thus 
results in a maximum allowable crystal temperature of 23.5°K. With the 
allowable ZiT = 18.5°K, solution of equation (23) results in a lead sink 
mass of 0^68 kg for the nominal mission described at the beginning of 
this section, which has a time interval of t = 2.57 x 10^ seconds. The 
required heat sink volume is 55 cm^. 

The heat sink mass and cold vacuum container will be 
highly polished to minimize the radiation heat load to the heat sink 
itself, and the cold mass will be supported within the container by 
means of thin stainless steel tension wires providing a negligible heat 
leak. 


(b) Collector Wall 

(1) Insulation Heat Leak - In order for the collector 
wall to function efficiently as a heat sink, it must be insulated from 
the surrounding vehicle structure which is nominally at 300 °K temperature. 
This is most conveniently accomplished by the use of multilayer aluminized 
mylar sheets, i.e., super insulation. The heat transfer through this 
insulation is by both conduction and radiation, and is commonly expressed 
as follows: 


Q - (5.67 X 10-12) . .^4) 


( 24 ) 


where 3 = 0.002 for optimum super insulation installation. 
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From the earlier discussion of capture coefficient 
criteria, the collector wall area, Ag, was related to the inlet aperture 
area, Aq, by the number of collisions required to reduce the incoming 
molecule temperature below a critical level, and to allow only 1 percent 
of the molecules to escape before being cryopumped. On this basis. 


A > 400A^. 
s o 


For the aforementioned crystal, i.e., A^, = 4 cm^, this 
criterion would imply Ag > 16 cm^. In this study a larger collector has 
been considered in order to reduce the escaping molecular fraction below 
1 percent and to arrive at a reasonable prototype design which groups all 
components, including the crystal heat sink mass, into a neatly packaged 
unit. A collector of 3.2 cm diameter by 8.9 cm length is therefore 
selected for the prototype design. This collector has a wall area of 
88.8 cm^. 


On the basis of this collector geometry, the collector 
wall heat load due to superinsulation heat leak is 


Q_ =8 milliwatts. 
^Ins 


(2) Gas Cooling Heat Load - The collector wall must 
cool the incoming mass flux to approximately 100°K before the gas 
molecules strike the crystal. The amount of energy which the wall must 
remove is essentially the kinetic energy of the incoming stream. Evalua- 
tion at 140 km altitude gives a conservative value of this energy removal 
requirement. 



= 84.5 milliwatts/cm^ 


(25) 


and 


Q , = 3.38 milliwatts, 
^cool 


(3) Collector Wall Heat Sink - The collector wall must 
serve not only as a heat sink but also structurally as a vacuum chamber 
containing the crystal sensor and the crystal heat sink. Aluminum is a 
more desirable metal in this application than the lead used in the crystal 
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heat sink, not only because of structural characteristics but also 
because of aluminum* s greater specific heat in the desired temperature 
range, i.e,, 70^K to 100°K. 

For the total collector wall heat load of 11,38 milli- 
watts, equation (23) indicates that 0,31 kg of aluminum is required to 
maintain the wall temperature below lOO'^K during the assumed mission. 

The required metal volume is 113 cm^, 

VI. SUMMARY DESCRIPTION OF RECCMIENDED SYSTEM 


On the basis of the greater sensitivity of the crystal sensor, the 
simpler coolant system, compared to those connected with the intermittent 
method, and the fact that an inlet shutter is not required, the continuous 
collection system is deemed the most promising instrument for prototype 
development. All major components of this system have been conservatively 
analyzed for an assumed flight, and reliable density measurements can be 
expected with an overall error (including telemetry) of only 6 percent. 

A summary description of the system follows. 

The densitometer consists of these major components: 

(1) A cylindrical orifice collector with a capture coef- 
ficient of 98 percent, 

(2) A cooled piezoelectric crystal which responds to mass 
addition from cryopumping with a frequency shift, 

(3) A measuring system consisting of a variable frequency 
oscillator and a beat frequency- to-voltage converter. 

(4) A refrigeration system consisting of two heat sinks, 
one of lead to cool the crystal to a temperature not 
higher than 23.5®K, and the other of aluminum which makes 
up the collector body and keeps itself at a temperature 
no higher than 100 ®K. 

Figure 9 shows these components, and Figure 7 depicts the measure- 
ment scheme. Aside from electric power and telemetry, the densitometer 
system can operate independently of the vehicle. Vacuum and temperature 
conditioning of the system would be achieved on the ground either in a 
fixed or a mobile laboratory with relatively simple equipment require- 
ments. Once conditioned and inserted into the vehicle, the system can 
withstand a hold of at least 12 hours without attention or umbilical 
connection. 
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Space, weight, power and telemetry requirements onboard a vehicle 
estimated as follows: 


(1) Space 

Dens itometer 
Electronics 


650 cm^ 
15,000 cm® 


(2) Weight 

Densitometer 

Electronics 


1.4 kg 

4.5 kg 


(3) Power - 30 watts at 28 volt DC. 

(4) Control Signals and Telemetry 

One control signal from onboard programmer to eject 
aperture cap. 

Four telemetry channels: 

a. 0-5 volt DC 

b. monitor variable oscillator setting 

c. lead heat sink temperature 

d. aluminum heat sink temperature. 
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Crystal Sensor Frequency Shift Measurement System 







APPENDIX A 


Capture Coefficient 


During the start of the collection process, the incident molecules 
will strike a cryopanel that will probably be covered with layers of 
absorbed gases. As the collection process continues, a condensate layer 
will form, and the incident molecules will impinge upon the condensate 
layer. In principle, it is necessary to estimate the sticking coef- 
ficient for the incident gas impinging on these two types of surfaces, 
but in practice it is expected that only the collisions of the molecules 
on the condensate will be important. Since the incident gas is air with 
a composition ccmiparable to the sea level condition, the important gases 
are nitrogen and oxygen. 

No experimental data have been found corresponding to a high 
velocity directed flow impinging on a cryopump. The available sticking 
coefficient data is limited to flows with S = 0 and temperature below 
400°K. Furthermore, investigators do not seem to be in agreement on 
the values of sticking coefficient. Wang et al. [8] found that the 
sticking coefficient increases as the thickness of the condensate layer 
increases. This increase ceased after a layer corresponding to about 
10"^ monolayers was deposited. This suggests that the collector plate 
sticking coefficient would vary with the thickness of the condensate 
layer; however, Buffham et al. [9] using what seems to be a comparable 
experimental arrangement found no variation of sticking coefficient with 
condensate layer thickness. There is also a lack of agreement concern- 
ing the level of the sticking coefficient. Dawson [10] presents the 
results of a large number of tests that determined the variation of 
sticking coefficient with both cryopump temperatures and tenq>eratures 
of the impinging gas. His sticking coefficients are significantly lower 
than the values reported by Buffham et al. [9], Bachler et al. [11] and 
Stickney and Dayton [12]. The experimental data are considered to be of 
limited value for two reasons: (1) the data do not cover the range of 

test conditions required, and (2) the data do not appear to be consistent 
among different investigations. In order to apply these data to the pre- 
sent analysis, theoretical methods must be used as a guide for extrapola- 
tion. 


The theoretical investigations of the sticking coefficient have 
generally been limited to considerations of the one -dimensional model 
consisting of a monatomic molecule colliding with a chain of elastically 
coupled lattice molecules. Since these correspond to normal collisions, 
the results might be suitable for the highly directed flow that exists 
in the present investigation. One objection to this approach is the use 
of the highly organized linear oscillator lattice for a condensate layer 
that appears to have many voids and defects, Roder [13]. Another 
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objection is the fact that the collision model is only correct for an 
infinite speed ratio, while the speed ratio for a 180 km orbit is about 
10. Although the random energy is small compared to the directed energy 
1/100), the random energy is comparable to the heat of vaporization 
of nitrogen and oxygen so the condensation process may be influenced by 
the finite speed ratio. 

Another theoretical approach has been considered by Buffham et al. 
[9]. They assumed that the gas arrived at the cryopump with a Maxwellian 
velocity distribution (S = 0) and all the molecules with an energy level 
below the critical energy level (E^) were captured. The sticking coef- 
ficient was then obtained from 


1 - 



(Ec/RT) 


(A-1) 


where 


T = temperature of the surrounding gas. 

Since no estimate was made for the value of the critical energy, it must 
be obtained from experiments. Figure A-1 is a plot of 1 - a vs 1/T for a 
range of values of the parameter E^. 

Included are data points for nitrogen. The data of Dawson [10] 
indicate that the critical energy level is about 1500 cal/mole which is 
comparable to the heat of vaporization. The other data include values 
in excess of 4000 cal/mole. Since most experiments were only accurate 
to 1 or 25 percent, there is considerable uncertainty in the critical 
energy to be associated with the sticking coefficients of nearly unity. 
This could be more readily resolved with data obtained at a higher gas 
temperature. 

From Figure A-1 it may be assumed that Dawson's data represent a 
lower limit of E^ ^ where E^ is the energy required to remove a 
molecule from the surface of the condensate. Since a prohibitive experi- 
mental accuracy would be required to set an upper limit on Eq using room 
temperature nitrogen, it is necessary to use a theoretical approach to 
estimate the upper limit. 

McCarroll and Ehrlich [14] used the one -dimensional lattice to 
estimate Eq/E ]3 (see Figure A-2). They found that, for equal molecular 
weights, the energy ratio Ec/Eb varied with the spring constant ratio, 
Kc/Kk* is the spring constant for the elastic coupling between the 

incident molecule and the lattice molecule, while is the spring con- 
stant between the lattice molecules. If the masses and the spring con- 
stants were the same, they found that E^/Eb ^ 25. 
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Only as approaches 0.2 does the energy ratio approach unity 

which is the value obtained from Dawson* s data. Since it is not known 
why this low value of Kc/K^ should exist, it appears that the energy 
ratio Ec/Eb can only be assumed to fall between 1 and 25. 

1.0 


1-0 0.1 


0.01 

0 .004 .008 .OlE .016 

1/T 

Figure A-1. Variations of Sticking Coefficient with 
Temperature 



““ " ' 






1 - r t T 1 

fc. 







Dawson __ 

Bachler et ah,- 
Stickney &. 

V V 






. ® 


\ 






. 1 

11^ 

i W 






1 

1 

1 1 







1 



Sl 





— ii 

/ay 


f 




\ 

\ 


Iv 












4( 

3< 

XX) 

yoc 

11 

17 

i 

50C 

XX) 

k 

E 

c 

- C 

a1 

/m< 

Die 




1 



1 

IX 










1 














1 


Cl 




\ 









1 


T~" 


\ 


\ 









\ 


"V 



L 


k 








\ 







5 








T 



L 


1 



S, 







J 

1 


\ 



V 









1 

[ 


n 



> 

\ 



\ 





aaaQ[Aa/>Qw'Q 


K. K 




Figure A-2. One-Dimensional Model of Surface Interaction 
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APPENDIX B 


Accommodation Coefficient 


Since the energy level of the incoming molecules is about two orders 
of magnitude greater than the heat of sublimation while the critical 
energy for condensation is about the same order of magnitude as the heat 
of sublimation, it is clear that the incoming molecules must be cooled 
by successive collisions with the cold wall before condensation can 
occur. The number of collisions required will depend upon the accommoda* 
tion coefficient. 


a = 


E. 

X 

E. 


X 


^off 


wall 


(B-1) 


The accommodation coefficient varies with such parameters as 

(1) gas species 

(2) wall species 

(3) gas temperature 

(4) wall temperatures 

(5) contamination of the wall surface. 

There are no experimental data on the accommodation coefficient 
for high velocity, highly directional air impinging on a cryogenic sur- 
face which is covered with a layer of condensate. 

Some limited data were obtained by Mayer et al. [15] using a cold 
plate normal to a supersonic stream. Since the flow was not free molec- 
ular relative to the plate, the temperature of the incident gas was lower 
than the stagnation temperature, and the data are not applicable to the 
present analysis. 

In general, molecular beam data are not suitable since they have been 
obtained at relatively low energy levels and with no condensate layer. 


Theoretical results often suggest the hard sphere approximation for 
high energy levels. 


a = 


(Ml + Ms) 2 




(B-2) 


where 


M = molecular weight. 

Unfortunately, this simple approximation is of questionable value 
if the molecular weight of the gas is equal to that of the solid, since 
it implies complete accommodation. The approach of Oman et al. [16] 
would undoubtedly be useful, but their results have not included the 
case where = Ms. 

Rogers [17] did include the case of Mx = Ms in his investigations 
using an analog computer, but since his **wall" consisted of a single 
linear oscillator, it is questionable whether the results are quantitative. 
Other investigations have found that the effects of the collision penetrate 
into the wall for several molecular diameters during the time required for 
the collision to be completed, McCarroll and Ehrlich [14] have included 
calculations of the accommodation coefficient for Mx = Ms over a range of 
Kc/Kk (see Figure B-1). 

Using the results of Appendix A, it is assumed that and 

Kc/Kk = 0.2; therefore, from Figure B-1, it appears that a = 0.77 should 
be cons erva tive . 



Figure B-1. Accommodation Coefficient Variation 
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APPENDIX C 


Sputtering 


When high velocity molecules or ions impinge on the surface of a 
target, they may cause some of the target material to be ejected or 
sputtered. This sputtering process is discussed quantitatively in terms 
of the yield which is defined as the ratio of the rate of sputtering of 
target molecules to the rate of incidence of energetic molecules or ions. 
The yield depends on the energy of the incident molecules, the binding 
energy of the target molecules, the angle of incidence and the atomic 
weights. 

It has generally been found that molecules or ions which are travel- 
ing at satellite velocities do not cause measurable sputtering on metals; 
thus, there should be no erosion of the basic condenser surface during 
the collection cycle. There is still the possibility that the incident 
molecules will cause sputtering of the condensate layer. This possibility 
is due to relatively low binding energy of the condensate layer when com- 
pared to that of the typical metal. 

Since no data were found on the self-sputtering of condensed nitrogen 
or oxygen, it was necessary to estimate the self-sputtering yield from the 
correlation of available data. 

From a comparison of the sputtering of metals by noble gas ions, it 
was found that the mass of the incident ion was not a dominant factor if 
the yield was low. Therefore, a crude correlation was made using the 
data of Stuart [18] for argon ions on various metals. The correlation is 
presented in Figure C-1, which is a plot of yield vs the ratio of the 
energy of the incident ion to the heat of sublimation (E^/E^). The plot 
shows that the data correlate within an order of magnitude. Assuming that 
nitrogen has a similar self-sputtering yield, and noting that E^/E^ ~ 
9/0.07 = 130, it appears that the yield will be about 1. If this is cor- 
rect, each incident molecule will sputter one molecule from the condensate 
layer. After the first collision, the energy level of the incident mole- 
cule is about (1 - CX) E^; therefore, if a = 0.77, the energy ratio 
is about 30. Since the sputtering associated with this energy ratio is 
small compared to the uncertainty in the initial sputtering yield, any 
sputtering after the initial collision can be neglected. 






APPENDIX D 


Refrigeration by Hydrogen Sublimation 


t 




A refrigeration system using the sublimation of solid hydrogen as 
the heat removal mechanism could be considered for cooling the crystal 
and the collector walls. This concept involves the relationship between 
the equilibrium vapor pressure of a substance and the substance tenq>era- 
ture. A heat load impressed upon the substance can be removed from the 
system under constant temperature conditions if the substance sublimes 
and the evolved vapor is removed from the system at a rate which main- 
tains a constant equilibrium vapor pressure. The heat removal capacity 
of such a technique is equivalent to the heat of sublimation of the sub- 
stance which is equal to the sum of the heat of fusion and the heat of 
vaporization at the particular temperature of interest. A survey of 
candidate substances for such a scheme indicates that only solid phase 
neon and hydrogen meet the primary requirement of solid phase temperature 
near 20®K at readily controllable vapor pressure, e.g., greater than 
1 torr (Figure D-1). The heat of sublimation of hydrogen (at 13.95°K) 
is 


/SB, - 125 cal/gm 


(D-1) 


P torr 



M . 



Figure D-1. Solid Phase Vapor Pressure 
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as compared to the heat of sublimation of neon (at 16®K) of 


^Ne “ cal/gm. 


(D-2) 


Therefore, on a weight basis the hydrogen system is definitely superior. 
On the other hand, the density of solid neon is much greater than solid 
hydrogen, and as a result a greater volume of solid hydrogen is required 
for a given heat load than would be required for solid neon. 


hydrogen (^^Hg) p = 10.6 cal/cm® (D-3) 

neon P“ 36.1 cal/cm^. (D-4) 


A secondary consideration is the vapor pressure level implied by the 
operating conditions for the two coolants. The vapor pressure of solid 
hydrogen at the triple point is 


Pv„ = 54 torr, (T = 13.95“K) (D-5) 

which is the highest possible operating temperature for this particular 
coolant. The vapor pressure of neon at the same temperature is 

= 0.1 torr, (T = 13.95‘*K) (D-6) 

and indicates that considerably better vacuum would be required for the 
use of neon as a coolant at this temperature. The crystal sensor, how- 
ever, can operate successfully as a mass collector at temperatures as 
high as 20®K, and since the triple point of neon is above this value 
(T’pp = 24.57 ®K), the neon could be operated at this higher temperature. 
In this case, the required vapor pressure could be 


P« = 28 torr. (T = 20®K) (D-7) 

*Ne 

A somewhat better vacuum condition is still required of the neon coolant 
system as compared to the hydrogen system. 



Hydrogen is the most promising substance for this sublimation scheme 
on the basis of the weight advantage, and because of cost considerations 
which include the availability at the launch site. In order to avoid a 
double coolant system, i.e., liquid nitrogen cooling for the collector 
walls and hydrogen sublimation for crystal cooling, the solid hydrogen 
volume can be sized to provide sufficient cooling for both functions. 

The heat of sublimation plus the crystal heat load implies a 
hydrogen sublimation rate which in all probability is not matched to 
the required gaseous heat removal capacity in the 70°K or 100°K tempera- 
ture range required of the collector walls. Consequently, the greater 
of the two heat loads will size the mass of solid hydrogen needed for 
a given flight time, and thus the coolant system will have excess capacity 
for the lesser heat load. This heat capacity matching poses problems in 
both situations; i.e., if the collector wall heat load requires the 
greatest hydrogen gas flow, a heat load beyond that supplied by the 
crystal must be imposed on the solid hydrogen to maintain the desired 
sublimation rate. On the other hand, if the crystal heat load is high, 
the gas flow to the collector may be more than required; this would 
result in a low collector wall temperature which could begin to act as 
a condenser, thereby invalidating the crystal sensor measurement. If 
such is the case, an artificial heat load must be imposed on the subliming 
gas before it reaches the collector section. 

From Section V, the sample gas cooling load plus the heat load due 
to superinsulation was calculated to be 11.38 milliwatts. Hydrogen gas 
provides 80 cal/gm cooling capacity between 70®K and 100®K; thus, the 
required hydrogen sublimation rate for collector cooling is 


collector 


. _ 11.38 milliwatt 
® ~ 80 cal/gm 


3.39 X 10"^ gm/sec. (D-8) 


The crystal heat load from Section V was 1.037 milliwatts with 
negligible heat loss due to insulation. The required sublimation rate 
for the crystal is thus 


crys tal 


. _ 1.037 milliwatt _ 
^ 125 cal/gm 


gm/sec. 


(D-9) 



The collector cooling therefore requires the greater hydrogen gas 
flow; however, there should be no problem in providing the desired subli- 
mation rate by proper design of the superinsulation section surrounding 
the solid hydrogen (Figure D-2). 


Mr> warmed to lOO^K 
/ ^ 
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Quartz Crystal Microbalance 


A. 




Figure D-2, Hydrogen Sublimation Coolant System 


The sublimed gas can be circulated through an intermediate layer 
of the superinsulation thus providing the proper temperature gradient in 
the super insulation for the desired sublimation rate, and at the same 
time, warming the hydrogen to 70°K before entering the collector coolant 
section. Tb^ required volume of solid hydrogen for a given flight time 
is 


V T . j tt (cm^) = 
solid H 2 


3.39 X 10 


-5 


= 3.94 X 10“^ t. 


(D-10) 


j 
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For the hypothetical mission considered in this report t = 3.57 x 10^ 
seconds and the required volume is 140.5 cm^. 

The maintenance of a constant vapor pressure over the solid hydrogen 
should pose no problems during the vehicle flight. The major heat loads 
are radiation loads and superinsulation losses and are invariant with 
time and vehicle altitude; thus, the sublimation rate will remain approxi- 
mately constant during the flight, and an orifice in the exhaust line can 
be sized to set the desired vapor pressure for the mass flow rate. 

Since the collector is evacuated before launch, the heat loads on 
the launch pad are essentially the same as the flight conditions. Main- 
taining the sub-atmospheric vapor pressure that is required if the hydrogen 
is to remain in a solid phase presents a problem, however, during any 
unattended hold periods, since the subliming hydrogen cannot be conven- 
iently exhausted by the atmosphere as it can in flight. Neither can the 
system simply be capped off after evacuation since a 10-liter volume would 
be required to hold the vapor pressure below 54 torr during a 12-hour hold 
period. 

An alternate approach would be a liquid fill and liquid phase boil- 
off to the atmosphere during the ground hold period followed by solidifica- 
tion at flight altitude as the liquid hydrogen is exposed to vacuum 
conditions. Adverse experience with hydrogen vent line plugging by 
freezing at high altitude would indicate that extreme care would need to 
be taken in the design and prototype development to avoid losing all the 
liquid before any solidification in the hydrogen container takes place. 

The most reliable hydrogen system would require continuous vacuum 
pumping during the ground hold period followed by orifice-controlled 
venting to vacuum while in flight. 
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